Anomalous Hall Effect and Skyrmion Number in Real- and Momentum-space 
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We study the anomalous Hall effect (AHE) for the double exchange model with the exchange cou- 
pling \Jh\ being smaller than the bandwidth \t\ for the purpose of clarifying the following unresolved 
and confusing issues: (i) the effect of the underlying lattice structure, (ii) the relation between AHE 
and the skyrmion number, (iii) the duality between real and momentum spaces, and (iv) the role 
of the disorder scatterings; which is more essential, an (Hall conductivity) or pu (Hall resistivity)? 
Starting from a generic expression for oh, we resolve all these issues and classify the regimes in the 
parameter space of Jht (t: elastic-scattering time), and A s (length scale of spin texture). There 
are two distinct mechanisms of AHE; one is characterized by the real-space skyrmion-number, and 
the other by momentum-space skyrmion-density at the Fermi level, which work in different regimes 
of the parameter space. 

PACS numbers: 72.15.Eb,75.47.-m,75.50.Pp,75.70.-i 



The anomalous Hall effect (AHE) is a phenomenon 
where the Hall resistivity has an additional contribution 
due to the spontaneous magnetization in ferromagnets. 
This anomalous contribution has been attributed to the 
spin-orbit interaction, and various mechanisms has been 
proposed H 0- Recently it has been recognized 
that the original expression by Karplus and Luttingcr 
0, i.e., the intrinsic contribution, has the geometrical 
meaning in terms of the Berry-phase curvature in mo- 
mentum space EJ 0- This is analogous to the the 
integer quantum Hall effect (fQHE) with the strong ex- 
ternal magnetic field |3, |9(. It was also proposed that 
AHE arises even without the spin-orbit interaction if the 
spin configuration is non-coplanar with finite spin chiral- 
ity, i.e., the solid angle subtended by the s pins where the 
electron hops successively [lol HH H3 . Uli. 1l4 | . Consider 
the double-exchange model 

H = ^2 t c\.Cr' - ^2 ' [ c r crc r] (1) 

(r,r') r 

where (r,r') runs the nearest neighbor sites, = 
(c^i , c^) ) is the annihilation (creation) operator at the 
site r, and S r is the classical spin localized at the site 
r. Assuming a strong Hund coupling |Jff|(>> \t\) be- 
tween the conduction electrons and the localized spins, 
the Berry phase of the localized spins acts as a fictitious 
magnetic field for the conduction electron 0, 0, fl7j |. 
Ye et al. assumed that this fictitious magnetic field has 
a uniform component due to the spin-orbit interaction 
in the presence of the uniform magnetization jfOj . How- 
ever there is a subtle issue concerning the definition of 
the skyrmion number when the spins are defined on the 
discrete points and/or the underlying lattice is relevant. 
This is related to the length scale with respect to the spin 
texture and/or the lattice structure. Furthermore, the ef- 



fect of the spin-orbit interaction can not be represented 
by the spatially uniform magnetic field; it induces the ef- 
fective "magnetic field", i.e., the Berry phase curvature, 
in momentum space. In real systems, the disorder is also 
relevant and often the following question arises: Which 
is more essential, the Hall conductivity an or the Hall 
resistivity pnl Therefore it is highly desirable to resolve 
all these issues in a unified fashion by clearly articulat- 
ing the connection between the AHE and the skyrmion 
number. In this paper, we study the AHE for the double 
exchange model eq. in the case of the small exchange 
coupling \Jh\ compared with the bandwidth \t\. 

It is found that the non-coplanar spin-configuration 
induces the AHE through the two distinct mechanisms 
originated by (Mi) the non-zero topological- windings of 
spin texture, and (Mn) the nontrivial structure of un- 
derlying lattice. Here the nontrivial structure means that 
the Wigner-Seitz unit cell contains multiple sites and dif- 
ferent kinds of plaquettes. Although the two mechanisms 
work simultaneously in generic cases, we can identify the 
two distinct limits where these mechanisms works exclu- 
sively: 

Ct: The unit cell contains a single site and the length 
scale A s of the spin texture is sufficiently longer than the 
underlying lattice constant a (A s ^> a). 
Cn- The unit cell contains multiple sites and the peri- 
odicity of the spin texture is the same as the underlying 
lattice (A s = a). 

The mechanism Mi is dominant and the AHE is char- 
acterized by the real-space skyrmion-number in the case 
Ci, while Mn is dominant and is characterized by the 
momentum-space skyrmion-density at the Fermi level in 
Cn. In the latter case, the AHE takes place even if the 
spin texture has no winding in real space and the total 
skyrmion-number is zero in momentum space. 

We also discuss the effect of disorder which causes the 
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following crossover. The Hall conductivity <jh is pro- 
portional to \Jfj\T 2 /\t\ for \J H \ < \t\(a/X s ) d < 1/r < 
\t\(a/X s ) in the case Ci, and for \Jh\ -C 1/r in the case 
Cii. Here r is the elastic-scattering time, and d is the 
dimensionality of a system. Hence the Hall resistivity 
Ph = -<th/(oq + ° 2 h) (o"o oc r, cr > |cth-|), does not 
depend on po — l/oo, where Co and po are diagonal 
conductivity and resistivity respectively. However, when 
1/r <C \Jh\i a H approaches its intrinsic value, and pn 
is proportional to p\ even in the case where | Jh | is suffi- 
ciently smaller than the band width. This suggests that 
the intrinsic (i.e., not due to impurity scatte ring ) mean- 
ing of the AHE can be observed through pn [la] even in 
the weak coupling regime. 

In order to consider the present issues with taking into 
account the effect of disorder, it is transparent to start 
with the Stfeda formula 1191. 
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x < Tr 



dEf F (E) 
dG 

J ll -^-(E)J v G s (E)+R.c. 



•(2) 



where V the system volume, J the current operator, 
G±(E) = [E-H- H imp ± G s = G+ - GL, H imp 

being the interaction with impurities, and (• • • ); mp repre- 
sents the ensemble average over impurity configuration. 
Here and hereafter we take the unit where c = h = 1. 

In this paper, we shall focus on only metallic states 
(|t|-7- > 1), and just replace G± by [E- H±i(2r)- 1 ]- 1 to 
bring in the effect of disorder, i.e. assuming the isotropic 
impurity potential and neglecting the localization effect. 
At zero temperature and using the symmetry Ov^^v — 
~^p,v\^vi the Hall conductivity is given by 



&H = 



+ 



V E B± - ■ l(^\J\lpa') X (l]} a ,\J\l/> a ) 

arctan(— 2r£ Q ) — arctan(— 2r£ a /) 

(£a-£a') 2 

2r(l+4r 2 £ a £ a Q 

(a-e ct 0[l + (2re Q ) 2 ][l + (2r^) 2 
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V 



J^e ± • [{ip a \J\t/} a ,) x (i]} a ,\J\il> a ) 



r 2 [f F (E a )-f F (E a ,)] 
l + [{E a -E a ,)T] 2 



(3) 



where £ Q = E a — uq (/xo : the chemical potential), 
fp(E) is Jf(E) with replacing the inverse temperature 
by 64r/(37r). The unit vector e± is normal to the plane 
determined by the Hall measurement. It is noted that 
cry depends on e±. The second formula in eq. J3J is an 
approximation for the first one by using physical impli- 
cation, and it is numerically confirmed that this approx- 
imation is almost exact with the accuracy less than 5% 



around peaks. In the clean limit, r — > oo, eq. J2J reduces 
to the formula obtained in Refs. [1,0. 

Now we shall consider the AHE in the weak-coupling 
regime by assuming a periodic configuration of localized 
spins. Every lattice-site is classified into sublattices I = 
A, B,C ■ ■ ■ . The localized spins satisfy S r< zi — Si. Then 
the coupling with the localized spins, H' , is expressed as 



H' 



where cit* is a Fourier transformation of the annihila- 
tion (creation) operator at I-sublattice, c^jj. The sys- 
tem has multiple bands even if we start with a single 
band at Jh = 0. The band separation is of the order of 
\t\(a/ \ s ) / (kp\ s ) d ~ 1 except for the spin splitting of the 
order of \Jh\i where kp is the Fermi momentum in the 
case Jh = and d is the dimensionality of the system. 
For simplicity, we shall consider the case kp ~ l/a here- 
after. When the parameter Jh satisfies \ Jh\ \t\(a/X s ) d 
and k is not near level crossings, approximated eigenval- 
ues and eigenfunctions are obtained by the conventional 
perturbation theory for spin-degenerate systems, 
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I=A,B,C,- 
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crcik 



(4) 



E n k± 
S n k 



E n k T -^-\S n k\ 



(5) 
(6) 



Si\u nkj \ 2 , 

I=A,B,C,— 

IV'nfei) = E E E C [n'k S ]lnk±]\^nL) + ' * ' » ( 7 ) 
n' s=± 1=0 

E U nkjXnk±^c\ k JO), (8) 

I=A,B,C,-<T=1,l 



where k s '][nks] * s ^ ne P er turbative coefficient of 
the Z-th order, u nk = (u nk ,A, u nk .B, u nkiC , ■ ■ ■ ) T is 
the orbital part of an eigenfunction in Jh = 0, and 
Xnk± = (Xnk±,-\,Xnk±,l) T is the spin part of an 
eigenfunction satisfying [S„k • &}Xnk± = ±\S n k\Xnk±- 
It is noted that S n k is regarded as the effective spin 
which is felt by the n-th pair of bands. Here, we 
use the word "the n-th pair of bands" for the bands 
with indices [nfc+] and [nk— ]. Below we shall consider 
the case Ci and Cn separately based on the above results. 

Case Ci 

In this case, the orbital part of an eigenfunction is given 
by u nk = u n = 1/Vflmb(l, e ib "- aB , e lb - ac , ■ ■ ■ ) T , where 
A^sub is the number of sublattices, b n is a reciprocal lattice 
vector, ai is a lattice vector between J-sublattice and A- 
sublattice. It is noted that u nk has no fc-dependence, and 
both S nk and Xnk± do not depend on the index [nk] in 
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this case. Using eqs. an is estimated as 

IstBZ / 



<JH = - 



VE nk X VE n , k 



+VE n , k x VE n „ k + VE n „ k x VE nk 
Jfj^R [S nn > ■ (S n ' n " x S n >> n )] 

^{E n k — E n 'k)(E n 'k — E n »k){E n i'k — E n k) 

y r 2 [F F {E nk )- r F {E n , k )] 
l + [{E nk -E n , k ) T ] 2 ' 



(9) 



Where S nn > = J2l U *nk,I S I U n'k,I = J2l U *n,I S I U n'J: and 

the intra-pair contribution (n — n') has been neglected 
because it is of the order of J H . The symbol J2' n n ' n" 
represents the summation with the condition that every 
index is different from each other. More precisely, S nn ' 
in eq. 10 is expressed as S[ nk ][ n > k ], and S nk in eq. © 



is Si 



nk] [nk] • 



However, in the case Ci, Sr. 



nk] [n' k] 



has 



fc-dependence. 

Changing the relative scales of the parameters \t\, \Jh\, 
1/X S and 1/r, er# shows the crossover. Here we shall 
consider the following two typical cases. 
Ci-A: | J H | < \t\(a/X s ) d < 1/r < |i|(a/A s ). This 
means that the length scale A s is shorter than the elastic- 
scattering length £ ~ \t\ar, i.e. A s -c £■ In this case 
\<th\ oc \J H \ 3 r 2 /\t\. 

Ci-B: 1/r <C | Jh\ <C \t\(a/X s ) d , i.e., the level separation 
is sufficiently larger than 1/r, where an takes its intrinsic 
value. 

In the former case Ci-A, especially in two dimensions, 
we can derive the topological meaning of the AHE by 
relating aji to the skyrmion number. Because of the 
condition \ Jh\ *C \t\(a/X s ) d -C 1/r, we can approximate 
the summations of band indices in eq. © by energy in- 
tegrals, and estimated them by residues E — no ± z/(2r) 
in the complex energy-plane. (See the first formula of 
eq. ©•) By the inverse Fourier transformation to real- 
space variables, we can obtain the expression equivalent 
to that given by Tatara and Kawamura 20] for the pe- 
riodic spin-configuration. They identified that \Jh\t is 
the small parameter for the perturbative expansion of 
<jh- Hence finite r is essential there. The Hall conduc- 
tivity was shown to be proportional to a sum of the spin 
chirality of any three localized spins with a geometrical 
weight in real space. 
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where x is the total chirality in real space, 
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I = r\ — r2, I' = T2 — r*3 and I" — r 3 — n. I(r) oc 
Efc e lfe r /[l + (2r£ fc ) 2 ] is an RKKY-type function which 
decay exponentially in the scale of I because of the 
complex part of the residues in energy integrals, and 
I'(r) — ^p-- Contribution from largely separated three 
spins with the scale of r decays rapidly as ~ e~ 3r / 2e , and 
the AHE is dominantly driven by chiralities of spins on 
small triangles. Using the expansion S ri = S r2 +i-VS r2 , 
S r3 = S r . 2 — I' ■ VS r2 , we obtain 



X 



S 3 A 



J d 2 r$ z (r), 



where S 



\S r \, N is the number of lattice sites, 
1 



$ z (r) 



AnS 3 



s ■ (y x s x v y s) 



is the skyrmion density in real space, and 



.4 



1 



dl 



dlTl' 2 l'(l)l'{l')I(l"). 



(12) 



(13) 



(14) 



Because of A s -C I, A is a dimensionless function of X s /a, 
when higher order terms in X s /£ is neglected. Thus the 
Hall conductivity is given by 

„2 73 ^2 c3 A AT 

" ' (15) 



er J 3 h t 2 



N 



where N s = J d 2 r<fr z (r) is the skyrmion number |2lll2 
Case Gn 

Here the unit cell contains multiple sites and the pe- 
riodicity of spin texture is the same as the underlying 
lattice. In this case, the band separation is of the order 
of \t\ except for the spin degeneracy and the degener- 
acy at the symmetric points. The orbital part of the 
eigenfunction u nk has non-trivial fc-dependence, and the 
same goes for S nk and Xnk- As we shall see below, the 
intra-pair contribution, which contains matrix element 
(ipnk T \J\if>nk±), is of the order of | J H \ 3 T 2 /\t\ or \ J H \/\t\ 
depending on | Jh\t "C 1 or 3> 1. The inter-pair contribu- 
tion, which contains the matrix element (ip n ks I J|VVfes') 
(n 7^ n'), is of the order of | J#| 3 /|i| 3 , and thus negligibly 
small compared with the intra-pair contribution. There- 
fore, the dominant contributions to <jh come from the 
intra-pair terms in sharp contrast to the case Ci. 

Considering that the n-th pair of bands is intersecting 
the Fermi level, the dominant contribution is given by 

- . IstBZ 
It \ s 

k 



[(lpnk+\J\lp,ik-) X (l/> nk -\J\lp. 



nk+) 



1-2 [fF^nk+) - fpiEnk-)] 



(16) 



1 + [{Enk- - Enk+)r] 

The perturbative eigenfunctions give the following matrix 
element of the current operator 

(lf>nk T \J\ll>nk±) = -^JH\S n k\xlk T [V fc e„fc ■ Cr] Xnfe±, 

(17) 
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where the higher order terms in Jh / 1 are neglected, and 

Cnfc = S n k/\S n k\. 

Using the above result and the completeness condition, 



t= 1.0, 6 = 71/3 



E 

*=+,- 



'X.nks 



1 X n fcs 



1 o 
1 



(18) 



we can estimate ur as follows, 



IstBZ 
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2V ^ 



~{E nk ) 



e i <I> 



n k 



(J H \S nk \) 3 T 2 
l + (J H \S nk \Tf 



where 



(19) 



(20) 



It is noted that $> nk may be regarded as the effective 
chirality of the n-th pair of bands in momentum space, 
because its integration is related to the solid angle of S nk - 
Especially in a two-dimensional system, e±_ ■ (47r) 
represent the skyrmion density in momentum space. This 
means that an is characterized by the momentum-space 
skyrmion-density at the Fermi level. When we change 
the parameter |</#|r, an shows the crossover, 



\Jh\ 



\uh\ oc 



\Jh\ 
1*1 



\Jh\t<. 1) 
|7ff|r» 1) 



(21) 



Finally, it is noted that, in contrast to the case Ci, the 
AHE takes place even if there is no winding of spin tex- 
ture in real space and the momentum-space skyrmion- 
number is also zero. 

In order to confirm the above consideration for the case 
Cn, we present the explicit results in Fig.^for the model 
cq. Jy on the kagome lattice with the spin texture, 



Sj_ 

s 



(sin 9 cos 4>i , sin 9 sin <pj , cos ( 



(22) 



where 4>j = (4nj — l)ic/6 (tia — 1, hb = 2, np = 3). The 
results are calculated by using the first line of eq. J3J), 
and clearly shows the crossover as predicted in eq. (|21|l . 

In conclusion, we have shown that the non-coplanar 
spin configuration induces the AHE by two distinct mech- 
anisms. The AHE is characterized by the real-space 
skyrmion-number when the underlying lattice structure 
is irrelevant and the spin texture is slowly varying as 
a <C A s <C L On the other hand, the AHE is charac- 
terized by the momentum-space skyrmion-density at the 
Fermi level when the underlying lattice structure is rele- 
vant and the periodicity of the spin texture is the same as 
the lattice. The Hall resistivity and conductivity become 
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FIG. 1: (a) r-dependence of pn and an, and (b) ^in- 
dependence of oh for the kagome-lattice model in the case 
where the average number of electrons per unit cell is three. 
It is shown in the panel (a) that pu and oh approach to 
constant values for \Jh\t ^ 1 and for \Jh\t S> 1 respec- 



tively. The panel (b) shows the crossover from \oh\ oc 
for \Jh\t < 1 to \oh\ oc \Jh\ for \Jh\t » I. 



Jh 



essential, i.e., independent of the elastic-scattering time, 
for | Jh\t <C 1 and for \Jh\t 1 respectively. 
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